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Abstract — The implications of causality are examined in detail 
for the permittivity, permeability, and refractive index of passive 
and active media. Particular attention is given to negative index 
media. 



I. Introduction 

The ability of fabricating artificial materials with sophisti- 
cated electromagnetic properties have generated large interest 
recently. By tailoring the electric permittivity e and the mag- 
netic permeability fi, one can obtain arbitrary metamaterial- 
induced coordinate transformations, leading to devices such 
as perfect lenses and electromagnetic cloaks [1|-|3|. 

At a single, positive frequency, a linear medium can have 
any e and p. If the imaginary parts of one or both of these 
parameters are negative, the medium is active; otherwise it 
may be passive. When considering global frequency-dependent 
parameters e(uj) and there are certain implications and 

limitations due to causality. We will analyze these connections 
in more detail, and summarize several results about passive and 
active media, beyond the usual Kramers-Kronig relations. We 
emphasize that only causality (and passitivity in the case of a 
passive medium) are taken into account. Thus the results are 
fundamental, and not related to practical fabrication problems. 

We restrict the attention to linear, isotropic, homogeneous, 
time-shift-invariant media without spatial dispersion. 

II. Passive media 

The permittivity of a passive medium satisfies the Kramers- 
Kronig relations HP] 
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'if e(o;) has singularities at the real axis, the Kramers-Kronig relations 
must be modified. For example, if the medium is conducting at zero frequency, 
e(oj) is singular at ui = 0. Then, the Kramers-Kronig relations are retained if 
we subtract the singularity, i.e., make the substitution e(u) — > e(ui) — ia/w 
in {TJ, where a > is the zero frequency conductivity Bj. These media have 
larger loss for the same Ree(cj). 



where V denotes the Cauchy principal value. Moreover, the 
permittivity is defined for negative frequencies by the symme- 
try relation 

e(-w)=e», (2) 

so that their inverse Fourier transforms are real. In addition to 
{T]i-([2| we have: 

Ime(w) > for co > 0, (3) 

which is the passitivity condition. The losses, as given by the 
imaginary parts of the permittivity, can be vanishingly small; 
however they are always present unless we are considering 
vacuum Q. 

In the following we will assume that the permeability p(ui) 
also satisfies relations ([TJi-Q, and <j3j> for a passive medium, as 
is usually assumed in the literature. It is to be noted, however, 
that this is not always exactly true. Indeed, using a Kramers- 
Kronig relation analogously to (jTbJ to calculate Re p(0) from 
Im p{uS), we straightforwardly find that Re /i(0) > 1 whenever 
Im/i(cj) > 0. This is not always correct |5), as demonstrated 
by the existence of passive, diamagnetic media. 

Eqs. ([T|l-(|3]l imply that e(uj) is analytic and zero-free in 
the upper half-plane Imw > p). Thus the refractive index 
n(u>) = W e{uj)^J ^{ui) can always be chosen as an analytic 
function there. With the additional choice that n(ui) +1 
as ui —> oo, n(ui) is determined uniquely, and it follows that 
(|TJ-([3j hold for the substitution e(w) -> n(u>) |6), (7]. 

We will now consider the implications of causality and 
passitivity. We will write the conditions in terms of the 
refractive index n(ui), although identical conditions apply to 
the permittivity and the permeability as well. 

A. Negative refraction is possible at a single frequency, with 
arbitrarily low maximum loss for all u>. 

The proof is by construction. Let the refractive index be 
written n(ui) = 1 + u(u>) + iv(ui), where 1 + u(lu) and v(u>) 
are the real and imaginary parts, respectively. Passitivity means 
that v(oj) > for oj > 0. S ince v(ui) = can be approached, 
we allow ourselves to put v(ui) = as well; adding a small, 
slowly varying function to v(lj) does not alter the argument 
below. Let v(ui) > vq > for ujq < u> < u>± — Alj, and 




Fig. 1. The loss function v(ui) has support below ui\. For uiq < to < 
ui — Acj, we assume v(uj) > vg. 



v(lo) = for u) > u)\, see Fig. [T| Estimating by the 

Kramers-Kronig relation (jTbJ, we find 
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Now we can choose any small vq and require the maximum 
of v(oo) to be less than, say 2vq- By having a sufficiently 
narrow transition region Aw, it(wi) gets less than — 1, which 
completes the proof. 

Note that the required function can be approached by 
superpositions of several, narrow Lorentzians, with resonance 
frequencies equally spaced in the interval (wq,Wi — Aw). 
For example, in the limit of continuous varying resonance 
frequencies from wq = to wi — Aw, we obtain 



it(w) + iv(u>) oc 
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where s = iy/uj 2 + iujT). The required refractive index 
function is obtained by choosing a sufficiently small T. 

The steep edge in the transition band may imply that the 
medium is difficult to realize. Given a small maximum loss, 
v(w) < n max for all w, one can prove that such steep edges 
are the only ways to obtain negative refraction for passive 
media that satisfy Kramers-Kronig relations in the usual form 
([TJ. Indeed, for any square integrable function w(w) > with 
limited steepness (\dv(uj)/duj\ < 1/Aw for some Aw), 



m(wi) > 



In- 
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(6) 



Here wi is an observation frequency. The inequality (|6]l is 
found by a similar argument as that of Q, considering the fact 
that the least possible u(wi) is obtained when w(w) = v maK 
for < w < wi — Su, and v(w) decreases linearly to zero 
above W! — Scu. Here Su> is a positive parameter. For media 
with singularities of e(w) or //(w) at real frequencies (e.g. an 
ideal plasma), |6]l does not apply. We can conclude: 



B. Passive, negative refraction implies either ( i) large loss be- 
low the working frequency, or (ii) exponentially steep variation 
immediately below the working frequency, or (Hi) singularities 
at real frequencies. 

The trade-off between requirements (i) and (ii) is quantified 
by (|6}. 

C. Negative refraction is possible in any finite bandwidth, with 
arbitrarily low loss there. 

This fact follows from a similar argument as that of the 
estimation Q; the trick is to choose a loss function that is 
sufficiently large below the interval of interest, and zero or 
sufficiently small in and above the interval. 

D. If a medium is lossless in a finite bandwidth [wi,Wa]» me 
minimum variation of n in this bandwidth is 

|Ren(wi) - 1|(2A - A 2 ) for rc(wi) < 1. (7) 

Here we have defined the normalized bandwidth 

A = (w 2 - wi)/w 2 . (8) 

For n(wi) > 1 there is no lower bound for the variation of 

n(w). 

In particular, the derivative dn/dw is bounded from below: 



dn I 2\n{u)) - l|/w for n(w) < 1, 
dw }0 for n(w) > 1. 



The proof is relatively simple, based on the Kramers-Kronig 
relations |7J. 

E. IfRen(uj) is constant and less than unity in [wi, W2], there 
is a frequency in [wi , 0J2] with loss 



w 2 — w 2 

Imn(w) > |Ren(wi) - l - 

ZW1W2 



|Ren(wi) - 1|A + 0(A 2 ) 



(10) 



Thus, a constant Re n(w) in a finite bandwidth, less than unity, 
is not possible unless the loss is sufficiently large. The proof 
is somewhat complicated j7), (8). 

F. Consider a passive, so-called perfect lens |7y. If the lens is 
used for frequencies in a finite bandwidth [w\,w^, the spatial 
resolution is 2nd/\ ln(A/2)|. 

Here we have quantified the worst resolution with respect to 
frequency in the interval, and optimized the medium subject 
to ([T])-(pJ. Thus, for a linear improvement of the resolution, 
the bandwidth must be shrinked exponentially. The proof is 
given in Q. 



III. Active media 

The permittivity of an active medium satisfies (T) and (2]), 
but not ([3]). Similar relations are valid for the permeability. De- 
termining the sign of the refractive index n(u)) = \J e(cj)/i(w) 
is not straightforward, as it is for passive media. For example, 
for a certain frequency it is no longer true to say that causality 
dictates a certain sign of n, based on e and fi at this frequency. 
One really has to go back to first principles: Causality means 
that no signal can go faster than the speed of light in vacuum. 
Restricting ourselves to media with no poles or odd-order 
zeros of e(uj)fi(u>) in the upper half-plane, it can be shown that 
causality implies that n(ui) must be identified as an analytic 
function in the upper half-plane, with asymptotic behavior 
n(ui) — » +1 as oj — > 00 (To), (TTJ. This is in accordance 
with the classical result for passive media [12|. 

In practice, the sign of the refractive index is conveniently 
determined using the following result: 

A. If e(uj) [i(u}) is continuous and zero-free for real frequencies 
(except possibly at uj — OJ, and analytic and zero-free in the 
upper half-plane, tlie refractive index is given by 

n(w) = V|eHIIM^)|exp[i(^H + ^(w))/2)], (11) 

where <ys e (w) + (fn(uj) is the complex argument of e(ui)/j,(ui), 
unwrapped such that it is continuous for ui 7^ and such that 
it tends to as w —> ±00 (13). 

The phase unwrapping procedure indicates that the sign of 
the refractive index is not uniquely determined from e and fi 
at a single frequency; knowledge of the global functions e(w) 
and n(ui) are required to determine the sign, even at a single 
frequency: 

B. Two media with identical e and identical [i at a certain 
frequency may have different n at this frequency. 

Consider two nonmagnetic materials with ei(ui) = 1 + Xg( w ) 
and €2(^0) = [1 + Xr{u)] 2 , respectively, where Xr.g{w) are 
Lorentzians in the form 
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With the parameters F g 
for material 1, and F r ■ 



F r,g^r,g 



(12) 



-0.00016, T„ = 0.005, andw„ 



1 



3.7, T r = 0.005, and uj r = 0.6 for 



material 2, these two materials have ei^i^g) — 1 — 0.030z. 
With the help of (TTJ, we find the refractive indices n\{u) g ) « 
1 - 0.015? and n 2 (uj g ) w -1 + 0.015i. 

As /i = 1 in the above examples, we obtain the following 
result: 

C. There exist right-handed negative index media. 

Also, there exist left-handed positive index media. These facts 
were first demonstrated in Refs. fl4| , (15). It is useful to 
examine the right-handed negative index media in more detail. 
Let a unit-step modulated cosine wave be incident to a half- 
spaced filled with such a medium. The time-domain field after 
some time is given in Fig. [2] The field plot for small z shows 
the "backward wave", with phase velocity and Poynting's 
vector in the — z-direction. This wave draws energy from 
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Fig. 2. The electric field in a right-handed negative index medium. The 
medium fills the half-space z > 0, while there is vacuum for z < 0. A unit- 
step modulated cosine is incident from vacuum, and the plot is a snapshot at 
time t = 30. For details see Ref. |10| . Animations are given in Ref. |16| . 



the medium and grows in the —z-direction. As always, the 
precursor, located close to z = ct, moves at a speed c in the 
+z-direction. 

D. Define the squared longitudinal wavenumber k 2 = 



e(ui)iJL(ui)u] 2 — k 2 — k 2 . If e(uj)^i(u})uj 2 — k 2 — k 2 has poles 
or odd-order zeros in the upper half-plane, k z does not have 
physical meaning for real frequencies. 

As a special case, if e(ui)fi(ui) has poles or odd-order zeros 
in the upper half-plane, the refractive index n(uj) does not 
have physical meaning for real frequencies. Although such 
media may be causal, they do not obey the conventional 
Kramers-Rronig relations for real frequencies; we must let 
the frequencies be complex with imaginary values larger than 
those of all nonanalytic points of n. 

The above result follows from the fact that the physical 
fields are the real time-domain fields. To obtain the time- 
domain fields from the frequency-domain fields, one performs 
an inverse Fourier or Laplace transform, that is, an integral 
above all nonanalytic points in the upper half-plane. Per- 
forming the integral along the real axis leads necessarily to 
unphysical results as the integral becomes dependent on the 
choice of branch cut locations flO) . 

When constructing active metamaterials, one should exam- 
ine whether e(uj)p,(uj) has poles or odd-order zeros in the 
upper half-plane. If there are such zeros or poles, the medium 
has so-called absolute unstabilities, which makes it useless for 
linear, small-signal applications. An absolute instability means 
that the fields blow up at any fixed point in space; the pulse is 
not convected away as in conventional gain media |T0"| , fTT) , 

(n), CD- 

It is to be noted that for oblique incidence, even conven- 
tional, weak gain media (such as Erbium-doped glasses) may 
have nonanalytic points of k z in the upper half -plane ]T3) . 
Thus for oblique incidence one must be particularly careful 
with the use or interpretation of k z . 
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Fig. 3. Refractive index for the left-handed medium with gain compen- 
sation. At w = 1.06 we have lossless negative refraction, n = —1 and 
dlmra 2 (a;)/du; = 0. 

E. A causal refractive index function n{ui) can approximate 
any square integrable function /(w) on a finite bandwidth. 

The approximation may be achieved with any precision^} in 
the Krein-Nudel'man case fl9) at the expense of the norm 
of n(u>) — 1 outside the bandwidth of interest. Of course, 
a large norm outside the relevant bandwidth may imply 
difficulties of realization. Nevertheless, we note that causality 
does not prohibit e.g. n(ui) « — 1 with any precision, even 
in a finite bandwidth. The possibility of approximating any 
desired behavior in a limited bandwidth may seem useless 
unless the resulting medium is free from absolute instabilities. 
Fortunately, there is a remedy in that log f(uj) can be approxi- 
mated by a function g(ui) using Krein-Nudel'man, and setting 
e(uj) — exp[g(oj)] |20| . 

As an example of gain compensation of the losses associated 
with a left-handed resonance, consider the causal medium 
e{uj) = 1 + x r (w) + X s (w) and fx{ui) = 1 + x r (w), where 
Xr,g(u) are the Lorentzians defined in ( p"2| ). Taking ui r = 1, 
uj g = 1.058, F r = 0.25, F g = -0.0034, T r = 0.005, and 
T g = 0.02, we find that n{ui) = -1 and dlmn 2 (w)/dcj = 
for uj — 1.06. This medium has Ime(cj) < in the bandwidth 
[1.05,1.07]; thus it is net active there. This does not imply 
that the system is unstable unless the medium is infinite or 
put in a resonator configuration; examples of such stable 
systems include fiber optic amplifiers. The absence of absolute 
instabilities of this medium is guaranteed by the fact that e(w) 
and fi(uj) do not contain zeros or poles in the upper half- 
plane. We note that the gain compensation due to x g (w) has 
completely removed the loss which would have been present 
in the absence of this gain. 

IV. Conclusion 

For passive media, we have first stated that causality does 
not prohibit negative refraction with arbitrarily low losses, 

2 In mathematical terms, a function h £ H 2 satisfying h(— uj*) = h*{ui) 
can approximate any function / £ L 2 (ui , UI2) , as precisely as desired in 
the corresponding metric. Here H 2 denotes the Hardy space of the upper 
half-plane, and < uj\ < u>2 < 00. See e.g. Ref. |19| . 



even in a finite bandwidth. Negative refraction implies either 
(i) large loss below the working frequency, or (ii) exponentially 
steep variation immediately below the working frequency, or 
(iii) singularities at real frequencies. However, for a fixed ui 
there is a lower bound for the loss Imu(w) if dRe n(uj) / 'duj = 
0, or lower bound for |dRen(aj)/dcj| if Imre(w) = 0. 

For active media, there are no fundamental limitations in a 
finite bandwidth. Nevertheless, there are some peculiar results 
related to the sign of the refractive index, and backward waves. 
Two media with identical e and [i at a certain frequency ui 
can have different n. The sign of n(oj) must be determined 
from the global frequency-dependence of e(w) and If 
active media are to be used in linear applications, it is cruical 
to eliminate the presence of absolute instabilities; otherwise 
the refractive index does not have physical meaning for real 
frequencies. 
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